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abstract 


Transmission  of  infrared  radiation  in  the  8-13  and  16.  5-21  fjm  spectral 
regions  from  3 km  altitude  through  the  atmosphere  to  outer  space  has  been  eval- 
uated. Analysis  of  tropical  radiosonde  atmospheric  profiles  of  temperature 
and  humidity  shows  that  fluctuations  in  the  total  amount  of  water  vapor  above 
3 km  cause  the  transmission  to  vary  by  amounts  from  10%  to  two  orders  of 
magnitude.  Fluctuations  in  the  effective  temperature  or  pressure  for  absorp- 
tion are  shown  to  be  small.  The  Curtis -Godson  approximation  is  analyzed  in 
detail  for  this  problem.  Also,  effective  parameters  describing  an  atmospheric 
slant  path  are  formulated  for  absorption  by  the  water  vapor  continuum. 

The  total  amount  of  water  vapor  above  3 km  varies  chiefly  as  the 
trade  wind  inversion  rises  above  or  falls  below  that  altitude;  it  correlates 
poorly  with  the  absolute  humidity  at  3 km. 

The  atmospheric  transmission  can  be  obtained  in  real  time  from  meas- 
urements of  atmospheric  emission.  Fluctuations  in  the  effective  emission  tem- 
perature, calculated  from  radiosonde  data,  are  relatively  small.  However, 
further  study  of  the  effect  of  cirrus  clouds  on  the  transmission-emission  re- 
lationship would  be  desirable. 
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1.  0 INTRODUCTION  AND  SUMMARY 


■» 


Quantitative  interpretation  of  infrared  measurements  through  the 
atmosphere  requires  accurate  specification  of  atmospheric  absorption. 
Convenient  yet  accurate  models  of  radiative  transfer  are  needed,  as  well  as 
knowledge  of  the  concentration  of  atmospheric  absorbers.  Among  the  latter, 
water  vapor  is  notable  for  its  spatial  and  temporal  variability. 

For  the  past  three  and  one-half  years  a multichannel  infrared  tele- 
scope  located  at  an  altitude  of  3.  0 km  on  the  summit  of  Mt.  Haleakala  on 
the  island  of  Maui  in  the  Hawaiian  Islands  has  been  used  to  study  extra- 
atmospheric  sources.  This  instrument  is  mounted  on  a 1.2  m telescope  and 
has  two  spectral  channels  which  are  particularly  useful  for  measurements  on 
objects  of  low  color  temperature:  No.  5 at  8.  3 to  13.4  jjm  and  No.  6 at  16.  5 
to  20.  8 fjm.  The  site  lies  above  most  aerosols  and  cumulus  clouds.  Never- 
theless, atmospheric  absorption  is  significant  in  Channel  6 at  all  times  and 
in  Channel  5 under  humid  conditions  or  at  low  elevation  angles. 

In  this  report  we  show  that  the  fluctuations  in  the  amount  of  water 
vapor  above  3 km  in  tropical  atmospheres  are  sufficient  to  make  reference 
to  typical  or  model  conditions  of  little  value.  We  explore  the  feasibility  of 
obtaining  real-time  information  on  atmospheric  transmission  from  sky 
emission  measurements.  This  procedure  is  favored  by  our  demonstration 

The  instrument  is  known  as  the  Advanced  Multi-Color  Tracker  for  AMOS 
(AMTA),  and  is  a facility  of  the  ARPA-Maui  Optical  Station  (AMOS),  which 
is  operated  under  the  direction  of  the  Avco  Everett  Research  Laboratory,  Inc. 
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that  tho  effective  temperature  of  the  atmosphere  for  molecular  emission 
is  independent  of  elevation  angle  and  has  small  temporal  fluctuations. 

However,  there  are  some  discrepancies  between  our  theoretical  results 
and  field  measurements  of  sky  brightness  which  deserve  further  investigation, 
with  attention  to  the  optical  effects  of  thin  cirrus  clouds. 

This  work  includes  careful  formulation  of  effective  parameters  for 
determination  of  the  absorption  or  emission  along  a slant  path  through  the 
atmosphere.  Local  radiosonde  atmospheric  profiles  are  utilized  to  display 
the  temporal  fluctuations  of  these  parameters. 

In  Section  2.  0 we  show  that  the  extinction  properties  of  an  atmospheric 
slant  path  for  molecular  band  absorption  can  be  specified  with  an  accuracy  of 
a few  percent  by  three  integral  parameters.  These  are  the  optical  depth  and 
the  effective  strength  and  width  of  spectral  lines.  The  latter  two  parameters 
are  weighted  averages  along  the  path.  Effective  parameters  for  the  water 
vapor  continuum  are  also  developed  as  weighted  averages. 

The  theory  of  emission  from  molecular  constituents  in  the  atmosphere 
is  presented  in  Section  3.  0,  with  emphasis  on  the  calculation  of  the  effective 
atmospheric  temperature. 

In  Section  4.  0 we  discuss  tropical  meteorology,  and  display  profiles 
of  water  vapor  concentration  and  temperature  obtained  from  radiosonde  data 
for  the  Hawaiian  Islands.  The  history  of  the  effective  atmospheric  parameters 
calculated  from  these  tapes  is  displayed  for  a period  of  one  year.  The  amount 
of  water  vapor  above  an  altitude  of  3 km  is  shown  to  fluctuate  by  a factor  of 
50,  primarily  because  of  variations  of  the  altitude  of  the  trade  wind  inversion 
or  because  of  passage  of  large-scale  cyclonic  disturbances  during  the  winter 


months.  The  amount  of  water  vapor  above  3 km  is  poorly  correlated  with 
the  local  water  vapor  concentration  at  3 km.  On  the  other  hand,  the  effective 
strengths  and  widths  of  absorbing  spectral  lines  show  very  little  temporal 
fluctuation.  Also,  the  effective  temperature  for  atmospheric  emission  varies 
relatively  little  over  the  course  of  a year.  Consequently,  all  atmospheric 
properties  which  determine  molecular  absorption  or  emission  above  3 km 
may  be  predicted  with  reasonable  accuracy  in  advance  except  for  the  amount 
of  water  vapor  overhead. 

Predictions  of  atmospheric  transmission  in  the  8-13  and  16.  5-21  4 rn 
spectral  intervals  are  presented  in  Section  5.  0,  with  elevation  angle  and 
amount  of  water  vapor  above  3 km  as  parameters.  The  effective  width  of 
each  spectral  line  is  calculated,  using  the  AFCRL  digital  compilation  of 
spectral  line  parameters.  The  transmission  averaged  over  several  lines 
is  then  obtained  from  appropriate  band  models.  This  scheme  yields  values 
of  transmission  by  water  vapor  which  are  somewhat  higher  than  laboratory 
measurements  at  Johns  Hopkins,  although  the  spectral  shapes  fit  well.  In 
the  16.  5-2  1 /im  spectral  region,  water  vapor  is  the  most  significant  absorber. 
For  low  elevations  and/or  conditions  of  high  humidity,  the  H^O  continuum 
(with  foreign  gas  effects  dominant)  absorbs  more  than  the  individual  lines. 
Transmission  predictions  based  upon  "average"  conditions  may  be  in  error 
by  as  much  as  an  order  of  magnitude.  The  8-13  Jim  window  is  much  more 
transparent.  Under  dry  conditions,  the  major  absorber  here  is  ozone,  whose 
effect  is  predictable  with  relative  confidence.  For  humid  conditions  the  H^O 
continuum  becomes  important.  As  the  water  vapor  concentration  increases, 
there  is  a transition  from  foreign  gas  effect  domination  to  self  "broadening". 
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The  overall  transmission  by  all  absorbing  species  is  presented  in 
Section  6.  The  dependence  of  transmission  on  elevation  angle  is  consistent 
with  results  from  the  single-parameter  LOWTRAN  3 code,  although  the 
latter  overestimates  the  transmission  in  the  16.  5-21  /im  region  by  as  much 
as  an  order  of  magnitude.  Some  discrepancies  are  found  between  our  trans- 
mission predictions  and  estimates  based  upon  atmospheric  emission  measure- 
ments. Further  analysis  of  the  relationship  between  these  two  should  include 
use  of  radiosonde  humidity  profiles  for  the  night  of  the  emission  measure- 
ments. Also,  the  emission,  absorption  and  scattering  by  cirrus  clouds 
should  be  addressed.  That  these  effects  may  be  significant  is  suggested  by 
comparison  of  predicted  and  measured  sky  radiance. 


i 
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2.  0 ATMOSPHERIC  ABSORPTION  - THEORY 


2.1  CURTIS -GODSON  APPROXIMATION 

We  will  show  that  to  calculate  the  spectral  absorption  along  a line 
of  sight  through  the  atmosphere  it  is  not  necessary  to  perform  a spatial 
integral  for  each  of  the  wavelengths  of  interest.  Rather,  the  atmospheric 
path  may  be  accurately  represented  by  only  three  integral  parameters 
which  relate  to  (1)  the  total  amount  of  absorbing  species,  (2)  the  effective 
average  temperature  along  the  path,  and  (3)  the  effective  average  pressure. 

This  approximation  applies  to  isolated  spectral  lines  or  to  a band  of 
overlapping  lines.  (For  continuum  absorption,  see  Section  2.  5.  ) Each  line 
may  be  conveniently  represented  by  its  equivalent  width,  defined  by 

W = J A(v  ) dv.  (1) 

line 

Here  A(v  ) is  the  spectral  absorptance  profile  of  the  line,  seen  at  high  spectral 
resolution.  It  is  shown  in  the  sketch  below  and  is  given  by 


A(v  ) = 1 - exp  - / k(v  ) ndjg  . 

path 


(2) 


The  absorption  coefficient,  k(v  ),  and  the  concentration  of  absorbing  species, 
n,  vary  along  the  path  with  distance  I 
For  a band  of  lines,  the 
spectral  absorption  averaged  over  a 
region  A v containing  many  lines  is 
of  interest: 


A = 1 - exp  - J kj(v  ) ndf  (3) 

path  i 

The  summation  is  over  all  lines  which  contribute  to  absorption  in  Av  . 

If  the  lines  are  weak  with  any  spectral  arrangement,  or  are  strong  and 
arranged  at  random  on  the  spectral  scale  (statistical  model),  the  average 
absorption  is  given  by 

X = 1 - exp  - "W/d,  (4) 

where  W and  d are  the  average  equivalent  width  and  average  spectral 
spacing  within  Av.  If  the  lines  are  strong,  of  the  same  strength,  and 
regularly  spaced  (Elsasser  model),  we  have 


A = 


(5) 


These  relations  are  plotted  in  Figure  1. 

Combining  Eqs.  (1)  and  (2)  for  a Lorentz  profile,  we  obtain 


W = 


/ 1 

line  ' 


1 - exp  - 


/ \ 

path 


(v 


S v du 

T ' 


dv . 


V 


(6) 


+ V 


(exact) 

Here  S = J k(v)  dv  is  the  integrated  absorption  coefficient  ("strength")  of 
the  line,  y is  its  half -width  at  half  maximum,  and  vQ  is  the  wavenumber 
at  its  center. 

The  Curtis -Godson  approximation  describes  the  absorbing  path  by 
an  equivalent  homogeneous  gas  sample  having  effective  line  strength  S, 
effective  line  half -width  y , and  optical  depth 


path 


/ n d|. 
path 


(7) 
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Figure  1. 


Average  Absorption,  Emissivity  and  Transmission  of  a 
Spectral  Band 


TRANSMISSION 


For  such  a sample  Eq.  (6)  becomes 


-f  ; s*;r  .i}*- 

lin.  V ^•V0)  + \ ' 


(C.G.)  (8) 


This  integral  has  been  evaluated  numerically.  The  result  may  be  written 


W = 2 w y f (x) 


(C.G.)  (9) 


where  f(x)  is  the  Ladenburg-Reiche  function  which  appears  in  Figure  2, 


T-  (K 

We  seek  S and  7 such  that  the  approximation  becomes  exact  in  the 
weak  (x<.2)  and  strong  (x>3)  line  limits.  For  weak  lines  the  exponent  in 
Eq.  (6)  is  very  small  relative  to  unity  for  all  wavenumbers,  and  the  exact 
equivalent  width  becomes 


•/(/  7 » >d;  1 ) 

line  ' path  v v ^ ^ 


l\L  ■ 


<v-vor 


■/ 

path 


S du 


(exact,  weak) 


For  weak  lines,  we  know  that  f(x)  ? X,  so  the  Curtis-Godson  approximation 
becomes,  from  Eqs.  (9)  and  (10), 


(C.G. , weak)  (12) 


Figure  2.  Curve  of  Growth  for  a Lorents  (Collision-Broadened)  Spectral 
e 
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For  the  approximation  to  be  exact  in  the  weak  line  limit,  conparison  of 
Eqs.  (11)  and  (12)  shows  that  the  following  relation  is  sufficient: 


S u = f S du, 
path 


or 


S = P*-th 


/ sn 


d£ 


f n d£ 
path 


(13) 


(14) 


(This  result  may  be  readily  obtained  for  any  line  shape.  ) Thus  the  Curtis- 
Godson  effective  line  strength  is  the  average  strength  along  the  path  with 
n,  the  concentration  of  the  absorbing  species,  as  a weighting  factor.  Since 
S depends  upon  temperature  but  not  pressure,  defines  an  effective  tem- 
perature for  the  absorbing  path. 

In  the  weak-line  limit  the  effective  line  width  y cancelled  out.  To 

determine  this  we  consider  the  strong-line  limit.  There,  the  exponential 

2 

in  Eq.  (6)  is  non-zero  only  for  values  of  (v  -v  ) much  larger  than  the 

2 

maximum  value  of  y (extreme  wings).  We  may  therefore  replace  y by 
y in  the  denominator  in  Eq.  (6): 


J " S 7 du 

W = J 1 1 - exp  - - - -P-fo- 
line  ' 


• ■ ■rTV 


dv 


(v-vor 


(15) 


(exact,  strong) 


Comparing  Eqs.  (8)  and  (15)  we  see  that  the  Curtis-Godson  approximation 
is  exact  in  the  strong -line  limit  if 


S y u 


f 

path 


S y du 


(16) 
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or 


7 


Y S n d( 


S n d( 


(17) 


The  Curtis -Godson  effective  line  width  ^ is  the  weighted  average  line  width, 
with  Sn  as  a weighting  function.  Because  7 is  proportional  to  pressure,  7 
might  be  considered  as  determining  an  effective  pressure  for  the  atmospheric 
column;  the  proportionality  of  7 to  the  inverse  square  root  of  temperature 
will  have  a minor  effect  in  ordinary  atmospheres. 

Along  a line  of  sight  at  an  angle  6 from  the  horizon  the  incremental 
optical  depth  in  any  stratified  atmosphere  may  be  written  du  = d d( 

= n esc  6 dh,  where  h is  the  height.  Thus  the  total  amount  of  absorber  is 


u = esc  9 J ndh  = _ q0o  esc  9 s esc  0 (18) 

path 

However,  in  the  C-G  parameters  £>  and  7 the  esc  Q factors  which  appear  in 
both  numerators  and  denominators  of  Eqs.  (14)  and  (17)  cancel,  so  that  S’ 
and  t*  are  independent  of  zenith  angle. 

Thus  the  C-G  parameters  may  all  be  expressed  in  terms  of  integrals 
with  respect  to  altitude.  The  weighting  functions  in  the  integrals  for  S and  y 
are  not  affected  by  the  elevation  angle  of  the  optical  path.  In  more  general 
terms,  the  weighting  is  not  affected  by  the  opacity  of  the  atmosphere.  Opacity 
effects  appear  in  the  use  of  the  C-G  parameters  to  calculate  atmospheric 
transmission,  as  in  Eqs.  (4),  (9)  and  (10). 
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2.2  ACCURACY  OF  C.  -G.  APPROXIMATION  IN  AN  EXPONENTIAL, 
MIXED,  ISOTHERMAL  MODEL  ATMOSPHERE 


The  Curtis -Godson  approximation  is  exact  in  the  weak  and  strong 
line  limits.  For  lines  of  intermediate  strength,  its  accuracy  may  be  checked 
for  certain  simple  model  atmospheres  for  which  exact  solutions  may  be 
readily  obtained. 

Let  us  consider  a model  atmosphere  described  as  follows: 


Total  concentration: 

Mole  fraction  of  absorber: 
Absorber  concentration: 
Temperature : 

Line  strength: 

Line  half -width: 


nt  = nt  (0)  e {//Ll 
X = const. 

n = X nt  = X nt(0)  e"*/Ll 


(19) 


:-f/Ll 


T = const. 

S = const. 

Y = Yo  nt  = \>  nt(0)  e' 

Here  the  argument  (0)  refers  to  any  reference  altitude  where  £ - 0. 
Lj  is  the  e-folding  distance  along  the  upward-looking  line  of  sight.  It  is  re- 
lated to  the  scale  height  Hj  by  Lj  = Hj  esc  9 where  6 is  the  elevation  angle. 
The  line  width  refers  to  Standard  Temperature  and  Pressure  (STP). 
Species  concentrations  are  to  be  measured  in  amagats,  i.  e.  , relative  to  the 
concentration  at  STP.  In  this  discussion  all  gases  are  considered  to  have 
the  same  broadening  power;  the  enhanced  broadening  power  of  H^O  is  negli- 
gible at  3 km  altitude  or  above. 

2.  2.  1 Curtis -Godson  Approximation 

The  effective  Curtis-Godson  parameters  for  a path  through  the  pre- 
sent model  atmosphere  to  outer  space  are  readily  calculated. 


oo  oo 

u = / n d£  = xnt(0)  j e‘*/Ll 

Jo  *'o 

= X nt  (0)  Lj 


dj£ 


(20) 


1 2- 


s = s 


V 


f v n df 
J 0 


= ^ont(0) 


= y ont(0)/2  = y ( 0 )/ Z . 


(21) 


(22) 


is 


Thus  in  the  Curtis -Godson  approximation  the  equivalent  line  width 
W = 2n  y f(x)  = 2 * f(2x1),  (C.G.  ) 


where 


x 


. su  _ , tlh 

' 2*  V ‘ 2tt  Yq 

SXL: 

X1  = 2jt  y ' 

o 


(23) 

(24) 


Although  the  parameters  u,  S and  y are  sufficient  for  this  analysis. 


it  is  instructive  to  introduce  two  additional  parameters:  the  effective  total 
density  ??t  and  effective  path  length  £,  defined  by  the  relations 

Y = Y0  srt  (25) 

and 

u = x nt  7- 

From  Eqs.  (20)  and  (22)  we  obtain 

TTt  = nt(0)/2  (26) 

and 


£ 


= 2 L, 


The  effective  total  concentration  (or  pressure)  is  one  -half  the  base -altitude 
value,  while  the  effective  path  length  is  twice  the  scale  length  Lj.  In  the 
amount  of  absorber,  u,  these  factors  cancel  (compare  Eqs,  (20  ) and  (25)). 
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It  should  be  observed  that  the  total  concentration  nt  does  not  appear 
in  Eq.  (24),  because  it  enters  into  both  the  numerator  and  denominator  of  the 
expression  for  x,  , 

For  a finite  path  through  this  model  atmosphere,  the  .generalization 
of  Eq.  (22)  leads  to  the  result  that  the  effective  pressure  for  a finite  path  is 
the  arithmetic  average  of  the  pressures  at  the  two  ends  of  the  path. 

2.  2.  2 Exact 

For  the  exponential,  mixed,  isothermal  model  atmosphere  Eq.  (6) 
become  s 


1 

7T 


2 -2</Ll  1 

s XV  0 [nt(0)p  e dg 

(v-Vo)2  + Yo2  lnt(0)lZ  «‘2£/Ll  - 


dv  . (27) 


Making  the  substitutions 

v - 1/ 

a = y (0)  ' P =2f/L1> 

and  using  Eq.  (22),  we  simplify  this  to 


W = y (0) 


exp  - 


da 


= 2Y(0) 


-Z  Y«» 


da 
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This  integral  has  been  evaluated  by  Strong  and  Plass:* 

l/2  Hx.  + !/2> 

W = 2 wl/£  v (0)  ' <Exact> 

The  gamma  function  T(x)  is  defined  for  all  positive  values  of  x. 
gral  and  half -integral  values  of  the  argument  we  have 
T(n)  = (n- 1 )!,  n = 1,  2,  3,  ..  . 

r(n  + 1/2)  = — - 5-  • • iZa~l)  r (1/2) 

_n 


(28) 

For  inte- 


r(l/2)  = 

The  asymptotes  of  the  gamma  function  are  such  that  Eqs.  (23)  and  (28)  agree 
for  small  and  large  x^: 


W — 27ry(0)  x^,  Xj  -*•  0 (Weak) 

W — 2y  (0)  (ttXj)1^2,  Xj  — ® (Strong) 


(29) 


The  approximate  and  exact  results  for  this  model  atmosphere  are 
compared  in  Figure  3.  The  match  is  excellent.  The  precentage  difference 
in  effective  line  width  W,  given  in  Figure  4,  peaks  at  6%  for  x^  = 0.  5 (x  = 1). 

The  6%  error  will  generally  be  an  upper  limit.  Water  vapor,  which 
is  of  prime  concern  to  us,  has  significant  concentration  over  a smaller 
pressure  range  than  the  uniformly  mixed  absorber  considered  above,  so 
the  Curtis -Godson  approximation  will  be  more  accurate  than  6%.  Also,  the 
lines  of  concern  to  us  are  rather  black  (x  >3,  Xj  > 6),  so  the  expected 
error  will  be  only  - 1%.  The  errbr  may  be  somewhat  higher  in  the  case 
of  absorption  by  ozone,  which  is  spread  over  a 20-km  thick  layer  --  twice 
a scale  height. 
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«7sie  Xj 


Figure  4.  Error  of  Curtis-Godson  Approximation 
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2.3  C.  -C,.  APPROXIMATION  APPLIED  TO  EXPONENTIALLY  MIXED. 
ISOTHERMAL  ATMOSPHERE 


The  model  described  above  in  Section  2.  2 applies  to  absorption  by 
species  which  are  vertically  mixed,  such  as  CO^  up  to  50  km  and  N2<D  and 
CO  up  to  10  km  altitude. 

The  mole  fraction  of  water  vapor  generally  falls  off  with  increasing 
height  in  the  troposphere,  where  the  variation  may  be  idealized  by  an 
exponential: 


"l/L, 

X = x(°)  e 

-ft/  L3 

n = X nt  = n(0)  e 

where 


(30) 


1 


Several  model  atmospheres  have  been  presented  in  the  AFCRL 
Handbook,  ^ ^ and  quoted  by  McClatchey.  The  scale  heights  H (corre- 
sponding to  scale  lengths  L for  a vertical  path)  for  these  models  are  given 
in  Table  1.  They  represent  the  altitude  increment  above  a reference 
altitude  hQ  in  which  the  quantity  shown  falls  by  a factor  of  1 /e.  In  these 
models  the  mixing  is  not  exactly  exponential,  as  shown  by  the  variation 
of  scale  height  with  altitude  for  the  tropical  model. 

For  isothermal  model  atmosphere  in  which  the  mole  fraction  of 
absorber  varies  exponentially  the  Curtis-Godson  parameters  are  readily 
obtained: 
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TABLE  1.  SCALE  HEIGHTS  IN  MODEL  ATMOSPHERES 
(AFCRL  Handbook,  1965) 


Reference 


Model  altitude 

hQ(km) 

U.  S.  Standard  0. 

Mid-latitude 

Winter  0. 

Mid-latitude 
Summer  0. 

Tropical  0. 


Total  density 

i 

9. 

8. 

9.  6 
9. 

9.  2 1 


H70  density 

2.  6 

2.  8 

2.  2 
2.4 

9 1.6 


1/e  Height  (km) 
^O  mole  fraction 

3 3.  6 

8 4.  1 

2.9 

8 3.  2 


3. 


(31) 


L1L2 

"<°>  ^3  = n(0)  l-tt; 


i 


S = S(const.  ) 

L.  L.+L- 

? = Y(o»  r-nq  ■ v(0) 


Also,  the  use  of  Eq.  (25)  yields 


nt  = 


UTT-  V°>  • 


If  we  let  — ■ a0 1 the  previous  results  for  uniform  mixing  (Eqs. 

(20)-(22)  and  (26))  are  obtained. 

The  exact  solution  for  this  model  atmosphere  could  be  obtained  by 
numerical  integrations,  which  we  have  not  attempted.  The  error  due  to 
the  use  of  the  Curtis-Godson  approximation  should  be  less  than  that  shown 
in  Figure  4 which  shows  a maximum  error  of  six  percent. 
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2.4  C.  -G  APPROXIMATION  APPLIED  TO  EXPONENTIALLY  MIXED 
ATMOSPHERE  WITH  LINEAR  TEMPERATURE  VARIATION 

Within  the  troposphere  the  temperature  decreases  with  height.  The 

variation  along  a line  of  sight  with  distance  f may  be  approximated  by  a 

linear  relation 


T = T 

o 

The  strength  of  an  atmospheric  line, 
temperature  T , is 

Q (T  ) /T  a"  r 

S(T)  = S(To,  (-f  ) [exp  - 


+ T'f  (32) 

relative  to  the  value  at  some  reference 


(33) 


o 

For  CO^,  n = 1;  for  H^O,  n = 1.5.  Also,  E"  is  the  lower  level  of  the  vi- 
brational transition.  At  tropospheric  temperatures  the  vibrational  partition 
functions  are  close  to  unity,  and  the  expression  exp  - is  very  small 
(for  v > 500  cm  *)  compared  with  unity.  Equation  (33)  may  therefore  be 
simplified  to 


S(T)  = S(To)(To/T)n  exp  - (^  . -j_)  (34) 

o 

We  now  evaluate  the  effect  of  this  temperature  variation  on  the 
Curtis -Godson  effective  line  strength  S in  an  atmosphere  in  which  the 
density  and  absorber  mole  fraction  decrease  exponentially  with  altitude 
and  the  temperature  decreases  linearly.  Substituting  Eqs.  (30)  and  (34) 
into  Eq.  (14)  yields 

/"  &)'  b - ^ (t  - ^)]b  - '/L’] « 
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The  upper  limit  is  the  distance  at  which  T — 0.  We  now  introduce  Eq.  (32) 
and  the  parameters 


A = hcE"/kT 
' o 

B = - T /L.  T' 
o'  3 


a = T/T 
' o 

to  obtain 

1 

S/Sc  = B exp(A-B)  J"  a”n  exp  - (—  - Ba)  da  (36) 

0 

This  expression  has  been  evaluated  numerically;  the  results  are  shown  in 
Table  2. 

The  relevant  values  of  the  parameters  may  be  obtained  for  the 
tropical  atmospheric  model.  There  we  find  that  T'  = -7  sin  Q (deg  K/km), 
vhere  Q is  the  elevation  angle.  Using  Tq  = 296°K  (the  reference  tempera - 
mre  for  the  AFCRL  compilation),  L3  = H3  esc  Q , and  H3  = 1.  6 km  we  have 

A = E"/206  cm-1 
B = 26 

Inspection  of  values  of  E"  for  major  HzO  lines  in  the  16-21  4 region  reveals 
-hat  most  of  them  lie  in  the  550-1100  cm-1  range,  so  that  A = 2.  7 to  5.  3 
covers  most  cases.  For  H,0  lines  Table  2 shows  that  S/S  = 0.  85  to  0.  95. 

L O 

This  result  indicates  that  the  temperature  dependence  of  the  line 
strengths  will  not  have  a large  effect  on  the  Curtis -Godson  parameter  S. 


-22- 
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TABLE  2 


n = 1 


CURTIS-  GODSON  EFFECTIVE  LINE  STRENGTH  S IN  A 
NONISOTHERMAL  ATMOSPHERE 


S/SQ  = Be(A'B)  J -L  e <Q  ^ da 
n a 


A 

B = 25 

35 

50 

70 

0.  1 

1. 040 

1.  030 

1. 020 

1. 016 

2.  5 

. 938 

. 959 

. 970 

. 982 

5 

. 853 

. 896 

. 924 

. 948 

7.  5 

. 784 

.841 

. 882 

.917 

0.  1 

1. 063 

1. 046 

1. 030 

I.  024 

2.  5 

. 957 

. 973 

. 980 

.989 

5 

. 869 

. 908 

. 933 

. 954 

7.  5 

. 796 

. 851 

. 890 

. 923 
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5 EFFECTIVE  PARAMETERS  FOR  THE  WATER  VAPOR  CONTINUUM 

The  Curtis-Godson  formalism  is  insufficient  to  describe  the  effect  of 
the  water  vapor  continuum.  Additional  parameters,  derived  in  the  present 
section,  are  required. 

The  atmospheric  transmission  by  the  water  vapor  continuum  at  any 
wavelength  may  be  written 


t 


exp  - 


/ kndl 

path 


(37) 


The  absorption  coefficient  depends  on  both  the  water  vapor  and  total  con- 


centrations: 


k = k (n  + a n.).  (38) 

o l 

The  coefficient  a varies  from  only  0.  001  in  the  8 to  12  ^m  window  to  about 
0.  02  at  20  /jxn.  ^ 

In  terms  of  the  absorber  mole  fraction  the  transmission  becomes 


t = exp  - / ko  (X  + ah  nt  dC-  (39) 

It  is  convenient  to  factor  t to  reflect  the  two  terms  of  k: 

4 = V 

where  the  effect  of  water  vapor  is  given  by 

tg  = exp  - / kon2df.  (40) 

and  the  effect  of  the  total  concentration  is  given  by 

tt  * exp  - / kQ  a n nt  d{  . (41) 

We  wish  to  obtain  effective  parameters  which  describe  an  equivalent 
uniform  gas  sample  having  the  same  transmission  as  the  atmospheric  path. 
First  we  treat  the  "seif"  factor  given  by  Eq.  (40).  If  the  temperature 
along  the  path  were  constant,  we  could  take  kQ  out  of  the  integral: 
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(42) 


ts  = exp  - kf 


dC 


. /n 

= exp  - kD  J"  n d u 

The  transmission  of  an  equivalent  sample  is 


t = exp  - k n u, 
s o 


(43) 


where  of  course  u = f n d C . Equating  these  expressions  for  tg  yields 

ff  u = J"  n d u 


or 

/n2  d* 

n - — 

J n df 


(44) 


Thus  the  effective  water  vapor  concentration  is  the  average  value  of  n 
along  the  path,  with  n also  appearing  as  a weighting  factor. 

Introduction  of  temperature  variation  along  the  path  leads  us  to  a 
third  effective  parameter,  kQ: 

t = exp  - k n u.  (45) 

s r o 

For  Eqs.  (40)  and  (45)  to  yield  the  same  result  we  must  have 


k n 
o 


u a 


n du 


or 


K »2  * 

f n-dl 


(46) 


For  this  parameter,  which  determines  an  effective  temperature  for  the 

2 

path,  the  weighting  factor  is  n . 

We  see  that  for  calculation  of  t any  atmospheric  path  may  be 
characterised  by  only  three  effective  parameters:  u,  n and  kQ.  The 
assumption  of  constant  atmospheric  temperature  was  merely  a temporary 
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expedient  to  guide  our  definition  of  n.  In  contrast  to  the  Curtis -Godson 
approach  to  a band  of  collision-broadened  lines,  no  approximation  is 
involved  here. 

An  effective  path  length  C could  be  used  as  an  alternate  to  n.  It 


is  naturally  defined  as 


7 = u/n 


[/ndt]' 
/ n2  d l 


In  terms  of  this  parameter,  Eq.  (45)  becomes 

ts  = exP  - kQ  u2/r.  (4 

To  obtain  effective  parameters  for  t^,  an  isothermal  atmosphere 
is  again  temporarily  invoked.  Equation  (41)  becomes 

tt  = exp  - kQ  q f nt  du 

(4 

= exp  - kQ  q nt  u 

The  effective  total  concentration  is  evidently  given  by 


"t  U = / nt 


/n  n d( 

* — -z — . (50) 

1 / n dj 

The  weighting  factor  is  the  same  as  for  if. 

A gas  sample  equivalent  to  a nonisothermal  atmosphere  has  trans- 


mission 


tt  ■ exp  - kQ  q «t  u 
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Comparing  Eq.  (41)  leads  us  to 


0T^)ntu=  J Yq  a nt 


du 


or 


(koa) 


yk  a n n df 
o t 

/ntndf 


Here  the  weighting  factor  is  the  product  nt  n. 


(52) 


These  parameters  may  be  evaulated  for  simple  model  atmospheres. 
If  the  temperature  and  mole  fraction  of  water  vapor  are  constant,  and  the 
total  gas  concentration  falls  off  exponentially,  as  given  by  Eq.  (19),  then 
evaluation  of  Eqs.  (44),  (47)  and  (50)  leads  to  the  results 

H * n (0)/2 

nt  = nt  (0)/2  (5  3) 


and 


« * 2Lj, 


in  agreement  with  Eq.  (26). 

If  the  mole  fraction  of  water  vapor  falls  off  exponentially,  as  in 
Eq.  (30),  we  obtain  readily 

n = n (0)/2 

T 

'1 


nt  = nt  (0)  Lj  + L3 


L.  + L< 


(54) 


x = = x <°)  — rr 

f = 2Lj. 

These  results  agree  with  Eq.  (31)  and  the  subsequent  equation. 
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3.0  ATMOSPHERIC  EMISSION- THEORY 


In  this  section  we  derive  relationships  between  absorption  and  emis- 
sion intensities,  with  special  emphasis  on  the  apparent  temperature  of  the 
atmosphere. 

3.  1 LINE  EMISSION 

The  atmospheric  radiance  (in  watts/ cm  - sr  - cm"  ) of  an  isolated 
spectral  line  may  be  written  compactly  as^ 


It  is  of  interest  to  evaluate  Eq.  (58)  under  limiting  conditions,  using 
the  Curtis -Godson  approximation.  In  the  weak-line  limit,  substitution  of 
Eq.  (12)  into  Eq.  (58)  yields 

f N°  (u)  S (u)  du 

N eff  = • (Weak  line)  (59) 

J ? (u)  du 

path 

o 

In  the  weak  line  limit,  Ng££  is  the  blackbody  radiance  averaged  over  dis- 
tance { along  the  path  with  S>  n as  a weighting  function.  Note  that  S (u) 
is  an  average  over  the  finite  path  of  optical  depth  u,  the  variable  of  integra- 
tion. 

In  the  strong  line  limit,  where  3?  » 1,  the  Ladenburg-Reiche  function 
has  the  asymptote  f(x)  ? (2x/jr)*/2.  From  Eqs.  (9)  and  (10)  we  then  have 

W (u)  = 2 (u)  y (u)  u]1/2.  (Strong  line)  (60) 

This  expression  might  be  differentiated  analytically  with  respect  to  u (see 
Ref.  5 for  such  an  approach),  but  the  result  is  complicated  because  of  the 
dependence  of  and  7 upon  u.  In  our  numerical  work  we  have  instead  cal- 
culated W vs  u using  Eq.  (60),  differentiated  the  result  with  respect  to  u, 
and  then  used  Eq.  (58)  directly. 

The  error  associated  with  the  use  of  the  Curtis -Godson  approximation 
in  determining  dW/du  may  be  greater  than  in  determining  W iteelf.  We  have 
not  evaluated  this  error.  It  may  be  that  the  criteria  for  "weak"  or  "strong" 
are  more  stringent  for  dW/du  than  for  W. 

In  the  strong  line  limit  one  may  be  concerned  that  the  integration 
indicated  in  Eq.  (58)  includes  small  values  of  u (u  < 2 it  y/S)  for  which 
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the  lines  are  weak  or  only  moderately  strong.  This  consideration  is  not 
a problem  if  (u)  is  essentially  constant  for  moderate  values  of  u. 

Then  a value  u^  » Ziry/S  may  be  defined  so  that  Eq.  (58)  becomes 

U1  u , 

/l  /*  end 

K S — / S ** 


= N 


/Uend 

K ^ du 


u. 


(61) 


where  uend  refers  to  the  end  of  the  path.  Here  W(Uj)  is  accurately  given 
by  the  strong  lire  limit,  even  though  the  details  of  the  function  dW/du  in 
the  range  0 < u < may  be  in  error  if  Eq.  (60)  is  used.  In  numerical 
cases  of  present  interest  we  have  determined  that  major  H20  lines  be- 
come strong  within  a fraction  of  a kilometer,  in  which  the  temperature  and 
hence  N°  are  indeed  constant. 

The  effective  temperature  of  the  atmosphere  (Teff)  for  line  emission 
may  be  calculated  from  N°^  by  inverting  the  Planck  function.  Equations  (58) 
to  (60)  show  that  T^^  may  be  obtained  from  atmospheric  temperature  profiles 
and  the  Curtis-Codson  parameters.  The  effective  blackbody  radiance  and 
temperature  of  the  atmosphere  for  line  radiation  are  relatively  insensitive 
to  the  elevation  angle  Q . For  example,  suppose  that  the  equivalent  line 
width  is  proportional  to  some  power  of  the  optical  depth  u:  W(u)  = Aun  = 

Au"  cscn0  . It  is  easy  to  show  that  the  integral  in  Eq.  (58)  is  also  proportional 
to  cscn0  , so  that  is  independent  of  9 . 

Actually  we  have  n = 1 for  small  values  of  u (weak-line  limit)  and 

n = \/Z  for  large  values  of  u (strong-line  limit).  But  in  practical  cases 

the  lines  are  weak  for  such  a small  range  of  u that  the  integral  in  Eq.  (61) 
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is  accurately  given  by  the  assumption  of  n = l/2  for  all  values  of  u,  as 

shown  by  Eq.  (61)  and  the  attendant  discussion. 

3.2  BAND  EMISSION 

Equation  (4)  tells  us  that  atmospheric  absorption  is  significant  only 
if  "W/d  (the  ratio  of  equivalent  line  width  to  spacing),  and  hence  the  effective 
amount  of  overlapping  of  spectral  lines  in  a band,  is  not  small.  Consequently, 
the  foregoing  analysis  of  emission  from  isolated  lines  needs  to  be  general- 
ized. The  following  is  based  in  part  on  Ref.  5. 

If  we  look  at  a band  with  spectral  resolution  that  is  fine  compared  to 
the  widths  of  the  lines,  the  spectral  radiance  of  the  gas  along  a path  through 
the  atmosphere  is 

N(  v ) = f N°  k(u)  t(u)  du,  (62) 

path 

where  the  k(u)  is  the  total  absorption  coefficient  due  to  all  lines  contributing 
at  wavenumber  v,  and  the  spectral  transmission  through  the  optical  depth 


u is 


Differentiation  yields 


t(u)  = exp  - j k(u')  du'. 
0 


= . t(u)  k(u) 


so  the  high  resolution  atmospheric  radiance  may  be  written  compactly 

N(  v ) = - J'  N°  dtd^  du.  (Band)  (63) 


path 
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averaged 


Of  more  interest  to  us  is  the  low  resolution  radiance  N 

v 

;r  some  spectral  interval  A v containing  many  spectral  lines; 


N = N(v)  = 


/ (/  N!  I du)  <wav 

A v 'path  / 


Inverting  the  order  of  integration,  and  noting  that  N°  is  constant  within 
A v , we  rewrite  this  expression  as 


Nv  = ' / Nv  (/  3^  dv/Av)  du 


N = 
v 


/ N°  <u)  dT  du' 
path 


where  t^  is  the  spectral  transmission  averaged  over  A 


ty  = y*  t dv/A  v • 


The  low  resolution  spectral  radiance  of  the  atmosphere  may  be  written, 


analogously  to  Eq.  (57),  as 


N = N ” € , 

v ell  v 


the  product  of  the  effective  spectral  blackbody  radiance  of  the  at- 
mosphere and  the  spectral  emissivity  of  the  atmospheric  path  e = 1 -t 
Equations  (64)  and  (65)  then  yield 


r _ dt 

/ N°  du 

rl .<  v d u 


(Band) 


This  expression  is  convenient  for  numerical  calculations.  N°„  is  the 

elf 
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blackbody  spectral  radiance  averaged  over  distance,  with  n dt  /du  as 

V 

weighting  function.  This  may  be  shown  more  explicitly: 

r n dt 

U n:  (u)  * du 

Th 


No  _ path 
eff  " 


(Band) 


(67) 


path 


du 


In  a spectral  band  described  by  the  statistical  model,  we  have 
t = exp  - W/d.  Substitution  into  Eq.  *66)  yields 


eu 


JL  N°  ^ H - W/J] du 

d (1  - exp  - W/d ) 


(68) 


If  the  effect  of  overlapping  is  negligible,  W/d  «1,  and  this  expression  re- 
duces to  Eq.  (58),  which  describes  the  nonoverlapping  line  case. 

For  a band  of  lines  we  have  shown  that  the  effective  blackbody  tem- 
perature of  the  atmosphere  may  be  calculated  from  temperature  profiles 
and  the  variation  of  atmospheric  transmission  along  the  line  of  sight. 

In  the  case  of  emission  from  isolated  lines,  we  saw  that  the  effective 
temperature  is  independent  of  elevation  angle  if  the  equivalent  line  width  W 
is  proportional  to  some  power  of  the  optical  depth  (u11).  The  angle  indepen- 
dence obtains  even  if  n is  dependent  upon  altitude,  provided  that  variation 
of  the  local  temperature  is  negligible  in  the  range  of  altitudes  in  which  n 
varies  (which  is  the  case  under  conditions  of  interest  here).  We  have  not 
succeeded  in  demonstrating  any  theorems  concerning  the  angle  independence 
of  effective  temperature  for  band  radiation. 
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4.  0 ABSORPTION  AND  EMISSION  PARAMETERS  FOR  A 
REAL  TROPICAL  ATMOSPHERE 


Infrared  transmission  through  the  atmosphere  is  significantly  affected 
by  temporal  fluctuations  in  atmospheric  conditions,  especially  water  vapor 
content.  The  shape  of  the  curve  of  water  vapor  mixing  ratio  vs  altitude 
often  departs  considerably  from  the  average.  Measurements  of  the  real 
atmosphere  must  be  considered. 

4.  1 TROPICAL  METEOROLOGY 

First  we  review  some  features  of  tropical  atmospheres,  with 
special  reference  to  the  Hawaiian  Islands  area.  Solar  heating  of  the  Earth's 
surface  near  the  equator  (intertropical  convergence  zone,  ITCZ)  produces  an 

upward  flow  of  air  which  then  travels  poleward  i ^ ^ t- n 1 1 it~ 

in  the  15-30°  latitude  belts,  which  include  the  Hawaiian  Islands.  Flow 
back  to  the  ITCZ  closes  the  loop.  The  downward  flow  velocity  is  typically 
60  m/day  at  the  700  mb  level,  falling  smoothly  to  zero  at  the  surface. 

At  the  higher  altitudes  this  flow  brings  down  dry  air  from  above. 

Also,  solar  heating  of  the  ocean  surface  below  produces  local  unstable  air 
masses  which  carry  moisture  upward  and  generate  cumulus  clouds.  Between 
these  masses  the  flow  is  downward.  The  combined  effect  of  these  flows  in 
the  lower  altitudes  is  a net  upward  transport  of  moisture  together  with  a 
net  downward  transport  of  air. 
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At  an  altitude  of  Z to  5 km,  the  upward  flow  of  moisture  and  the 
downward  flow  of  dry  air  meet,  producing  a sharp  step  in  the  profile  of 
absolute  humidity,  as  seen  in  Figures  5 and  6.  These  steps  are  almost 
always  at  the  same  altitude  as  temperature  minima  or  inversions,  shown 
in  Figure  7.  The  humidity  step  is  generally  stable  in  the  stratified  atmos- 
phere found  in  the  summer  (Figures  5 and  6),  when  the  meteorology  in  the 
Hawaiian  area  is  dominated  by  the  ENE  trade  wind;  hence  this  phenomenon 
is  termed  a trade-wind  inversion.  This  inversion  usually  marks  the  tops 
of  cumulus  clouds.  The  step  in  the  moisture  profile  is  preserved  because 
moisture  which  rises  up  into  the  inversion  layer  generally  condenses, 

(7) 

sometimes  in  the  form  of  "inversion"  stratus  or  stratocumulus  clouds. 
(Occasional  penetration  is  marked  by  "chimney  clouds".  ) Thus  the  upper 
dry  layer  stays  dry.  As  we  follow  the  trade  wind  from  the  Northeast 
Pacific  to  the  Hawaiian  Islands,  upward  transport  of  moisture  causes  the 
thickness  of  the  moist  layer  to  increase,  although  the  air  flow  is  slightly 
downward. 

This  trade-wind  dominated  pattern  also  occurs  at  times  in  the 
winter.  At  other  times  large  cyclonic  disturbances  dominate,  giving  rise 
to  much  greater  vertical  mixing.  Figure  8 thus  shows  humidity  steps 
in  some  profiles  and  widespread  mixing  or  patchiness  in  the  troposphere 
in  others  during  the  first  half  of  January. 
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TMTIAL  DATE=7/2/70 


FINAL  MTE=7/16/70 
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w*  ABSOLUTE  HUMIDITY  (GM/M*) 


Figure  5.  Humidity  Profiles  Derived  from  Rawinsonde  Data  for  Hilo, 
Hawaii  at  Midnight  GMT  (July) 


ABSOLUTE  HUMIDITY  (OM/M*) 


Figure  6.  Humidity  Profile*  for  Hilo  at  Noon  GMT  (July) 
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4.2  ATMOSPHERIC  RADIOSONDE  DATA 

Figures  5 to  8 illustrate  radiosonde  (balloon)  measurements  made 
twice  daily  at  0000  and  1200  hours  GMT  at  more  than  400  weather  stations 
throughout  the  world  and  available  on  magnetic  tapes  from  the  NOAA 
National  Climatic  Center,  Asheville,  North  Carolina.  The  tapes  give 
relative  humidity,  temperature,  and  wind  direction  and  speed  at  about  20 
altitudes.  We  have  analyzed  data  taken  at  Hilo  on  Hawaii  Island,  160  km 
SE  of  AMOS  Station  on  Maui. 

Some  caution  should  be  observed  in  using  the  tapes.  The  National 

(8) 

Weather  Service  hygristors'  ' have  reduced  sensitivity  at  relative  humidities 
less  than  20%: 

Relative  Humidity  (%) 

Actual  NWS  Hygristor  Reading 


> 20 

Actual  + 2% 

20 

20 

15 

18 

10 

16 

5 

14 

Above  the  trade-wind  inversion  altitude,  the  atmosphere  may  be  somewhat 
drier  than  indicated.  During  the  last  few  years  the  NWS  arbitrarily  replaced 
all  relative  humidity  readings  below  20%  by  the  default  value  of  19%.  The 
data  reported  here  is  from  earlier  years,  when  this  wrong -way  "correction" 
was  not  used. 
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4.  3 CLRTIS -GODSON  PARAMETERS  FOR  ABSORPTION  IN  HAWAIIAN 
ISLANDS  ATMOSPHERES 


Looking  at  Figures  5,  6 and  8 we  see  that  the  integrated  amount  of 
water  vapor, 


> 0 km 


- / 

3 km 


nHzO  dh* 


above  3 km  altitude  varies  over  a wide  range  of  values.  In  summer,  the 
variation  arises  as  the  trade-wind  inversion  rises  above  or  falls  below  the 
3 km  level.  In  winter,  storms  intermittently  raise  water  vapor  into  the 
upper  troposphere  which  is  otherwise  quite  dry. 

The  history  of  the  amount  of  water  vapor,  u,  above  3 km  at  Hilo 
is  shown  in  Figures  9 and  10  which  show  data  for  2 p.m.  and  2 a.  m.  local 
standard  time,  respectively.  The  extreme  values  differ  by  a factor  of  50. 
The  high  values  should  be  accurate  within  a few  percent,  but  the  low  values 
may  be  somewhat  high  due  to  hygristor  insensitivity. 

The  local  water  vapor  density  at  3 km  altitude,  p,  correlates  poorly 
with  the  total  amount  of  water  vapor  above  3 km,  as  shown  in  Figures  1 1 
and  12.  If  it  is  dry  locally,  it  is  dry  above;  but  if  it  is  humid  locally  there 
is  wide  scatter  in  the  total  water  above,  because  the  thickness  of  the  humid 
layer  above  3 km  varies  widely. 

The  other  two  Curtis -Godson  parameters,  the  effective  line  strength 
S defined  by  Eq.  (14)  and  the  effective  line  width  y defined  by  Eq.  (17)  , 
have  been  evaluated  using  Eq.  (34)  and  the  expression 

1/2 


V = Y0(P/P0)  (Tq/T)  *'  • 

Here  p is  the  total  atmospheric  pressure.  Self-broadening  may  be 

neglected  above  3 km.  Figures  13  and  14  show  S/SQ  and  V /y  Q for  the 

Hilo  radiosonde  data  above  3 km.  The  subscripts  ( ) refers  to  STP. 

o 


(69) 
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JULV  AUG.  SEPT  OCT.  NOV.  DEC.  JAN.  FEB . HAR.  APR.  HAY  JUNE 


Figure  9.  Amount  of  Water  Vapor  Above  3 km  for  Hilo  at  Midnight 
GMT  (1970-71) 
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Figure  1 1.  Total  Amount  of  Water  Vapor  Above  3 km  vs  Absolute 
Humidity  at  3 km  for  Hilo  (Summer) 
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V-  , 


WATER  VAPOR  ABOVE  3 KM, 
U (GM/CM»«2) 


Figure  12.  Total  Amount  of  Water  Vapor  Above  3 km  vs  Absolute 
Humidity  at  3 km  for  Hilo  (Winter) 
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m 


Figure  14.  History  of  Effective  Line  Strength  (S)  and  Width  (y  ) 
Above  3 km  for  Hilo  at  Noon  GMT  (1970-71) 


These  calculations  are  for  a lower  energy  level  E"  = 800  cm  *,  which 
is  typical  of  the  strong  lines  found  in  the  16-21  pm  window.  Varying 

E"  from  550  to  1100  cm"*  changes  S/Sq  and  y/y  only  by  1-2%. 

In  contrast  to  the  fluctuating  u,  the  parameters  S and  y are  seen 
to  have  very  small  diurnal  or  seasonal  variations  for  this  tropic*!  location 
because  the  amount  of  water  vapor  appears  in  both  the  numerator  and  de- 
nominator of  Eqs.  (14)  and  (17).  This  means  that  we  may  make  atmospheric 
absorption  calculations  assuming  constant  values  for  the  effective  tempera- 
ture and  pressure  of  the  atmosphere  above  3 km  at  Hilo.  (In  other  caleu- 
lations  for  temperate  atmospheres  S showed  a +_  10%  seasonal  change.  ) 

We  conclude  that  the  temporal  variations  in  the  atmospheric  column 
of  interest  may  be  described  by  changes  in  a single  parameter,  u,  as  far 
as  absorption  by  bands  of  water  vapor  lines  is  concerned.  For  the  water 
vapor  continuum  it  is  possible  that  additional  parameters  such  as  n or  nt> 
defined  by  Eqs.  (44)  and  (50),  might  be  necessary  for  complete  description 
of  local  meteorological  conditions. 
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4.  4 RELEVANCE  OF  HILO  DATA  TO  MT.  HALEAKALA 

Located  on  the  northeast  shore  of  the  Island  of  Hawaii,  Hilo  is 
generally  on  the  windward  side  of  Mauna  Loa  and  Mauna  Kea  (elevations 
4.  2 km).  (The  wind  at  1.  5 km  comes  only  rarely  from  more  than  + 45° 
away  from  the  east  - Ref.  9,  Figure  63.  ) Thus,  the  Hilo  rawinsonde  pro- 
files represent  the  interisland  trade  wind  atmosphere,  relatively  free 
from  orographic  influences. 

Radiative  heating  and  cooling  of  the  surface  of  the  mountains  (if  not 
hindered  by  cloud  cover)  generate  local  upslope  winds  by  day  and  downslope 
winds  by  night.  Mendonca^1^  has  studied  these  winds  at  Mauna  Loa  Ob- 
servatory, located  at  3.  4 km  altitude  on  the  windward  side  of  Mauna  Loa. 

He  found  the  nocturnal  downslope  wind,  which  lasts  until  about  0730  local 
time,  had  a thickness  of  only  about  50  m.  The  thickness  of  the  upslope 
wind  grew  to  some  500  m by  local  noon  - the  result  of  unstable  convective 
mixing  and  greater  fetch.  The  atmosphere  over  the  open  ocean  can  be  modi- 
fied close  to  a mountain  in  the  Hawaiian  Islands,  but  the  region  affected  will 
be  very  thin,  especially  at  the  top  of  a mountain  at  night. 

These  thermally  generated  local  mountain  winds  are  to  be  distin- 
guished from  the  mechanical  disturbance  of  the  trade-wind  flow  due  to  the 
mountains.  Because  their  slopes  are  very  steep,  the  trade  wind  will  gen- 
erally flow  around  Mauna  Loa  and  Mt.  Haleakala  more  than  over  them  (Ref.  6, 
pp.  105-6).  The  step  in  humidity  profile  has  been  found  at  the  surface  of 
Mauna  Loa  at  an  altitude  close  to  that  in  the  rawinsonde  data.  ^ ^ 
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4.  5 EMISSION  PARAMETERS  FOR  HAWAIIAN  ISLANDS  ATMOSPHERES 


In  Section  3.  1 we  showed  that  the  radiance  of  the  atmosphere  from 
an  isolated  spectral  line  ) can  be  obtained  as  the  product  of  the  effec- 

tive blackbody  radiance  (N^)  and  the  equivalent  line  width  (W).  The  last 
variable  depends  upon  the  absorption  parameters  u,  S and  y,  whose  varia- 
tions we  already  have  displayed.  Here,  we  consider  the  variability  of 


This  parameter  has  been  calculated  for  water  vapor  emission  using 
the  measured  atmospheric  profiles  3 km  and  more  above  Hilo,  Hawaii.  Nu- 
merical studies,  presented  below,  show  that  the  significant  emission  comes 
from  strong  lines,  so  we  use  Eqs.  (58)  and  (60);  for  comparison  we 

also  evaluate  Eq.  (59)  based  upon  the  weak  line  assumption.  These  results 
have  little  dependence  on  the  parameters  of  the  spectroscopic  lines.  The 
calculation  involves  determining  only  u,  ? and  y’,  and  then  W,  as  functions 
of  altitude,  and  we  have  already  seen  that  “S  and  7 are  affected  very  weakly 
by  the  choice  of  E",  the  lower  energy  level  of  the  spectral  line. 

Figure  15  shows  the  history  of  N°^,  the  effective  spectral  blackbody 
radiance  of  the  atmosphere,  at  1200  hours  GMT.  Inverting  the  Planck  func- 
tion yields  the  effective  atmospheric  temperature  histories  presented  in 
Figure  16.  In  each  figure  the  strong  line  assumption  gives  the  upper  curve, 
while  the  weak  line  assumption  gives  the  lower.  For  the  latter  we  "see" 
up  to  higher  altitudes  where  the  atmosphere  is  colder.  The  scales  of  these 
figures  are  quite  expanded;  a small  deviation  from  strong  line  conditions 
would  not  introduce  a significant  error. 

Referring  to  the  upper  curve  of  Figure  15,  we  see  that  N°  has  a 

eff 

seasonal  (summer  to  winter)  change  of  10%,  and  short-term  fluctuations 
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Figure  15.  Effective  Blackbody  Spectral  Radiance  for  Atmospheric 
Line  Emission  at  3 km  Above  Hilo 
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of  + 5%.  This  is  much  smaller  than  the  expected  fluctuation  of  W (propor- 
tional to  the  square  root  of  the  highly  variable  u)  or  N^n  . This  point  will 
be  elaborated  below. 

Similar  plots  for  0000  hrs  GMT  (2  pm  local  standard  time)  show 
values  of  Tfiff  generally  only  1-2°C  higher  than  for  1200  hrs  GMT.  The 
day-to-night  temperature  change  at  about  4 km  altitude  — the  significant 
level  — is  much  less  than  at  sea  level. 

Overlapping  of  H^O  lines  was  neglected  in  the  foregoing.  The  quan- 
tity W/d,  presented  in  Eq.  (4),  is  a measure  of  the  degree  of  overlapping. 

In  typical  numerical  cases  of  interest,  we  have  found  W/d  ~ 1 /2,  so  we  will 
now  consider  this  effect. 

The  spectral  atmospheric  radiance  of  a band  of  o/erlapping  lines  is 
given  by  Eqs.  (65)  and  (66).  These  cannot  be  evaluated  with  the  broad  gen- 
erality we  applied  to  line  radiation.  The  function  t^(u)  must  be  specified. 
Numerical  calculations,  described  below,  show  that  over  the  entire  range 
of  amounts  of  water  vapor  (u)  to  be  encountered  for  reasonable  zenith  angles 
the  transmission  averaged  over  the  17-21  ji  band  may  be  represented  as 

where  a 

sidering  that  the  absorption  is  dominated  by  strong  lines. 

Using  Eqs.  (66)  and  (70)  we  obtain  the  history  of  the  effect  blackbody 
radiance  of  the  band  shown  in  Figure  17,  and  the  effective  atmospheric 

temperature  in  Figure  18.  The  results  are  difficult  to  distinguish  from  the 
strong-line  curves  of  Figures  15  and  16.  The  effect  of  overlapping  has  a mini- 
mal effect  on  N°„  and  T ... 

eff  eff 


t (u)  = exp  - a u ' , (70 

1 /z 

= 0.  8 cm/gm  . This  type  of  dependence  is  to  be  expected  co>'- 
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Figure  18  Effective  Blackbody  Temperature  for  Atmospheric  band 
Radiation,  as  seen  from  3 km  Above  Hilo 


-56- 


The  zenith  angle  of  the  optical  path  does  not  affect  the  effective  black- 


body  radiance  for  isolated  lines;  if  W is  proportional  to  some  power  of  u, 
the  numerator  and  denominator  of  Eq.  (58)  change  with  angle  by  equal  fac- 
tors. Since  we  have  seen  that  overlapping  of  lines  affects  N°  very  little 

eff  1 

it  should  be  independent  of  zenith  angle  also  for  a band  of  lines.  We  have  not 
performed  detailed  numerical  tests  of  this  conclusion. 

The  spectral  radiance  of  atmospheric  H^O  lines,  averaged  over 
the  16-21  4m  region,  has  been  calculated  using  Eq,  (63),  The  result,  pre- 
sented in  Figure  19,  shows  considerable  day-to-day  fluctuations  - much  larger 
than  those  of  in  Figure  17.  This  means  that  the  atmospheric  H-»Q  band 

radiance  is  influenced  much  more  by  variations  in  the  amount  of  water  vapor 
overhead  than  by  variations  in  the  effective  temperature  of  the  water  vapor. 
Consequently,  if  we  can  perform  real-time  measurements  of  the  atmospheric 
emission  by  H^O  lines,  the  data  could  be  used  to  determine  e , by  assuming 
a constant  and  thus  the  water  vapor  line  transmission  t ^ = 1 - c . 

This  scheme  can  succeed  even  if  in  addition  to  H^O  line  emission 
we  have  radiation  from  the  continuum  or  from  other  sources,  because 

Eqs . (62)  - (66)  can  be  rewritten  in  terms  of  the  distance  ! instead  of  the 
optical  depth  u.  However,  cirrus  clouds  or  aerosols  which  do  not  obey  the 
relation  t = 1 - because  of  reflection  or  scattering  could  reduce  the 
accuracy  of  this  approach. 
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Figure  19.  Spectral  Band  Radiance  of  Atmospheric  Water  Vapor 
Above  3 km  at  Hilo 


-58- 


5.0  ABSORPTION  PREDICTIONS 


In  this  section  we  present  calculations  of  the  atmospheric  trans- 
mission along  lines  of  sight  extending  from  an  altitude  of  3 km  to  outer 
space,  with  elevation  angle  and  amount  of  water  vapor  overhead  as  var- 
iable parameters. 

5.1  SPECTROSCOPIC  PARAMETERS 

In  addition  to  the  Curtis -Godson  parameters  u,  §/Sq  and  y/yQ, 

obtained  in  Section  4 from  radiosonde  measurements  of  the  Hawaiian 

Islands  atmosphere,  we  use  the  AFCRL  compendium  of  spectroscopic 

line  parameters.^^  From  this  magnetic  tape  we  read  the  line  strength 

(S  ),  width  (v  ) and  wavenumber  for  each  spectral  line  of  the  species  of 
o o 

interest.  Transmission  by  bands  of  lines  was  calculated  using  Eqs.  (4), 

(9)  and  (10)  and  Figure  2. 

For  the  continuum,  we  use  the  data  of  Roberts,  Selby  and 

Biberman^  and  Burch,  Gryvnak  and  Gates.  ^ The  latter  indicates 
the  temperature  dependence  of  the  absorption  coefficient  kQ.  Equations 
(45),  (49)  and  (54)  were  used  to  evaluate  the  water  vapor  continuum  trans- 
mission. This  means  that  we  assumed  an  exponentially  mixed  atmospheric 
model  for  the  continuum;  the  l/e  fall -off  distance  for  H^O  mole  fraction  was 
estimated  from  Eq.  (31)  and  actual  Hilo  data  (Figures  11  and  12).  We  found 
H^  = 1.8  km,  which  is  close  to  the  value  (1.6  km)  for  the  tropical  model 
atmosphere  (see  Table  1). 

(12) 

The  AFCRL  spectroscopic  line  parameters'  ' for  H~0  have  been 
tested  against  measurements  at  John  Hopkins^ **  1 5 ^ of  absorption  along 
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6 passes  through  a 100  ft  long  tube  containing  water  vapor  and  nitrogen. 
Figures  20  and  21  compare  our  transmission  calculations  with  the  labora- 
tory data.  The  agreement  is  quite  good  in  the  20  - 25  (jm  case,  shown  in 
Figure  20.  Over  this  spectral  range  the  average  transmission  was  meas- 
ured to  be  59%,  while  the  calculation  yielded  57%. 

In  the  15  -20  pm  spectral  region  the  calculated  transmission  is 
generally  (25  - 50%)  higher  than  measured,  as  seen  in  Figure  21.  The 
lines  fit  fairly  well,  except  that  the  resolution  chosen  for  the  computed 
spectrum  (1  cm  * half  width  at  half  maximum)  is  a bit  high.  We  infer  that 
the  water  vapor  continuum  is  stronger  here  than  we  calculate  from  the  data 
of  Refs.  4 and  13. 

The  mole  fraction  of  water  vapor  is  0.  0027  (vs  a = . 02  to  . 04) 
in  the  case  shown  in  Figure  20,  while  % =0.04  (vs  a = .007  to  .028) 
for  that  shown  in  Figure  21.  Thus  the  absorption  coefficient  is  propor- 
tional to  the  total  concentration  in  the  first  case,  and  proportional  to  the 
H^O  concentration  in  the  second. 
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Figure  21.  Transmission  by  Water  Vapor  Over  600  ft  Path  (15.  2 - 20  ji) 
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5.  2 TRANSMISSION  IN  lb-21  gm  WINDOW 


Figure  22  shows  the  calculated  transmission  for  a vertical  line  of 
sight  originating  at  the  altitude  of  AMOS.  The  amount  of  water  vapor  over- 
head is  1 gm/cm^.  Figures  9 and  10  indicate  that  this  humid  condition  is 
encountered  occasionally  in  the  Hawaiian  Islands. 

The  major  absorber  is  seen  to  be  the  H^O  continuum  (t  = 30%).  The 
structure  in  the  spectrum  is  due  to  H^O  lines,  whose  transmission  averaged 
over  this  region  is  40%.  CO^  also  absorbs  in  the  16.  5-17.  5 pm  region. 

Ozone  absorption  is  negligible  at  these  wavelengths. 

The  transmission  by  the  individual  absorbing  components  has  been 

averaged  over  the  16.  5-20.  8 pm  region  and  displayed  in  Figures  23  to  26  for 

1 /2 

seven  cases.  For  the  water  vapor  lines  we  have  very  nearly  -lnt  o:  u ' , 
as  shown  in  Figure  23.  This  behavior  is  typical  of  strong  line  absorption. 
Indeed,  our  line-by-line  results  show  that  strong  lines  (x  i 3)  contribute 
80%  of  the  average  equivalent  line  width,  W,  even  for  the  smallest  value  of 
u considered. 

Individual  cases  for  the  water  vapor  continuum  calculations  are 
shown  as  data  points  in  Figure  23.  In  these  cases  the  total  gas  concentration, 
nt>  is  constant,  but  x.  nf  and  u = are  variable.  The  mole  fraction 
of  water  is  generally  small  compared  with  a,  the  total-concentration  coef- 
ficient appearing  in  Eqa,  (38)  and  (39).  We  thus  would  expect  Eq.  (49)  to  det- 
ermine the  transmission.  Accordingly  -lnt  should  scale  linearly  with  u. 

Our  individual  cases  fit  this  expectation  closely;  there  is  a little  scatter  in 
Figure  23  because  the  condition  x « a i«  not  perfectly  fulfilled  in  every 
case. 


-63- 


-In  (AVERAGE  TRANSMISSION) 


AMOUNT  OF  WATER  VAPOR  ALONG  LINE  OF  SIGHT,  u ( gm/cm2) 

GTatt 


Figure  23.  Transmission  by  Lines  and  Continuum  in  16-21  ^ Window  vs  u 
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Maps  of  water  vapor  transmission  versus  elevation  angle  and  amount 
of  water  vapor  overhead  appear  in  Figures  24  and  25.  Since  the  transmission 
depends  only  on  u,  the  constant  transmission  lines  shown  have  the  form 

u = u^  esc  6 - const. 

The  COj  transmission  depends  only  on  elevation  angle,  as  shown  in 

Figure  26.  Most  of  the  absorption  is  by  strong  lines  (x  > 10),  as  indicated 

1 /2 

by  the  - 1 nt  oc  (esc  8)  ' dependence  of  the  results. 

Figures  23  to  26  give  the  transmission  by  separate  atmospheric  con- 
stituents, each  averaged  over  the  16.  5-20.  8 jim  band.  The  average  trans- 
mission by  the  mixture  of  components  is  accurately  given  by  calculating 
the  product  of  the  individual  transmissions  at  each  wavelength  and  then  av- 
eraging over  the  spectral  band.  However,  the  overall  transmission  may  be 
obtained  to  an  accuracy  of  a few  percent  from  the  product  of  the  average  trans- 
missions. This  would  not  be  valid  if  the  positions  of  lines  of  different  species 
were  correlated  with  each  other. 

In  summary,  we  find  that  the  transmission  of  the  16-21  ^m  channel 
is  poor  (<  60%)  for  viewing  extra-atmospheric  sources  from  AMOS,  unless 
we  look  at  high  elevation  angles  under  unusually  dry  conditions. 
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Figure  26.  Transmission  by  CO^  Lines  vs  Elevation  Angle  in  16-21  jim  Band 
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5.  3 TRANSMISSION  IN  THE  8-13  gm  WINDOW 

The  8-13  jim  region  is  much  more  favorable  for  infrared  observations 
hrough  the  atmosphere.  Water  vapor  lines  are  relatively  transparent,  as 
seen  in  Figures  27  and  28.  The  significant  absorbing  lines  are  strong,  as 
shown  by  the  -lnt  oc  u ^ dependence  in  Figure  27,  but  they  are  widely 
separated  so  that  the  average  transmission  is  high.  For  small  amounts  of 
water  vapor.  Figure  27  shows  some  deviation  from  the  square-root  relation- 
ship as  the  lines  become  weaker. 

Our  results  for  individual  cases  calculated  for  the  water  vapor  con- 
tinuum are  shown  as  data  points  in  Figures  29  to  30.  The  mole  fraction  of 
water  vapor,  i®  indicated  for  each  case.  This  quantity  is  to  be  compared 
with  a,  the  total-gas  concentration  coefficient,  which  ra ages  from  - 0.001 
to  0.004  in  this  spectral  region.^’  1 ^ For  small  u (and  small  \)  we  have 
generally  \ « a,  so  the  transmission  is  approximately  given  by  Eq.  (49) 
and  -lnt  is  reasonably  proportional  to  u as  seen  in  Figure  29.  For  large 
u (and  large  \),  we  have  generally  \ » a,  so  the  transmission  is  approxi- 
mately given  by  Eq.  (48)  and  -lnt  is  proportional  to  u2/f  as  seen  in  Fig- 
ure 30. 

These  two  regimes  are  displayed  on  the  map  of  continuum  transmi  s ion 
(Figure  31).  For  small  lines  of  constant  transmission  are  given,  as  in 
•rigure  25,  by  the  relation  u = u^  esc  9 - const.  For  large  lines  of  con- 
stant transmission  are  given  by  u2 /f  - const.  Since  F,  defined  by  Eq.  (47), 

is  proportional  to  esc  9,  constant  transmission  corresponds  to  the  relation 

2 

u^  esc  9 - const. 

in  the  \ » a regime.  A more  detailed  analysis  would  produce  curvatu-  s to 
the  lines  in  Figure  31,  with  a smooth  transition  between  regimes. 
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Figure  27.  Transmission  by  H^O  Lines  in  8-13  (jm  Band  vs  u 
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Figure  28.  Map  of  Transmission  by  H.,0  Lines  in  8-13  jim  Band 
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Figure  29.  Transmission  by  H^O  Continuum  in  8-13  jim  Band  vs  u 


3- 


(AVERAGE  TRANSMISSION) 


-In  (AVERAGE  TRANSMISSION) 


.0002  0 


AVERAGE  TRANSMISSION 


CSC  (ELEVATION  ANGLE) 


AMOUNT  OF  WATER  VAPOR  OVERHEAD,  uA  (gm/sec2) 

07811 


Figure  31.  Map  of  Transmission  by 


H^O  Continuum  in  the  8-13  jjjm  Band 
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In  general,  absorption  by  water  vapor  in  the  8-13  pm  window  is  small 
under  dry  conditions,  but  is  significant  and  varies  rapidly  with  amounts  of 
water  vapor  under  humid  conditions. 

The  effects  of  CC^  and  are  presented  in  Figure  32.  Neither  is 
a serious  absorber  in  the  8-13  pm  window.  The  CO^  lines  are  of  interme- 
diate strength,  which  results  in  a -lnt  vs  esc  9 curve  which  is  intermediate 
between  a linear  and  square -root  relationship. 

The  significantly  absorbing  ozone  lines  are  strong,  and  at  individual 
wavelengths  we  find  that  -lnt  is  proportional  to  the  square  root  of  u or 
esc  9.  When  the  transmission  is  averaged  over  the  8.  3-13.4  pm  range,  a 
weaker  dependence  results,  as  seen  in  Figure  32. 
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Figure  32.  Transmission  by  CO^  and  in  the  8-13  /jim  Band 
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6.  0 CONCLUSIONS 


In  this  section  we  will  present  our  transmission  predictions  for  the 
Hawaiian  atmosphere  above  3 km  in  summary  form,  and  compare  them 
with  the  alternate  results. 

In  the  previous  section,  the  transmission  was  given  for  individual 
atmospheric  constituents  — H^O,  CO^  and  O^  lines,  as  well  as  the  H^O 
continuum.  For  any  amount  of  water  vapor  and  zenith  angle,  the  overall 
transmission  may  be  obtained  by  multiplication. 

Typical  results  for  our  detailed  calculations  of  the  overall  trans- 
mission are  shown  by  the  solid  lines  in  Figures  33  and  34.  The  lower  solid 
curves  are  for  = 0.  12  5 gm/cm2  of  water  vapor  overhead,  a "dry" 
condition  occasionally  encountered  at  night  (see  Figure  10).  The  upper  solid 
curves  are  for  n i = 1 gm/cm2  of  water  vapor  overhead.  This  humid 
condition  is  sometimes  found  day  or  night,  as  seen  in  Figures  9 and  10. 

The  slopes  of  these  curves  range  from  n = 0.  56  to  0.  74,  where 
n is  the  exponent  in  the  expression 

t = exp  — a (esc  0 )n.  (71) 

Figures  33  and  34  also  show  the  average  transmission  in  channels 
5 and  6 calculated  using  the  LOWTRAN  3 code16,  17  for  the  AFCRL 
tropical  model  atmosphere  (dash-dot  lines).  As  shown  on  Figures  9 and  10, 
this  model  entails  \x L = 0.  8 gm/cm  of  water  vapor  above  3 km  altitude. 

The  LOWTRAN  transmission  predictions,  relative  to  our  more  detailed 
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Figure  33.  Predicted  Transmission  from  3 km  to  Outer  Space  vs  Elevation 
Angle.  Present  results  for  large  (1  gm/cm2)  and  small 
(0.  125  gm/cm2)  amounts  of  water  vapor  overhead  are  compared 
with  LOWTRAN  calculations  and  predictions  based  upon  a 
measurement  of  atmospheric  radiance. 
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Figure  34.  Predicted  Transmission  from  3 km  to  Outer  Space  vs  Elevation 
Angle.  Present  results  for  large  (1  gm/cmZ)  and  small 
(0.  125  gm/cm2)  amounts  of  water  vapor  overhead  are  compared 
with  LOWTRAN  calculations  and  predictions  based  upon  a 
measurement  of  atmospheric  radiance. 
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calculations,  are  low  by  about  3 to  10%  in  the  8-13  /jm  region  and  high  by  fac- 
tors of  2 to  16  in  the  16.  5-2  lpm  region.  The  error  is  most  serious  in  the 
longer  wavelength  channel  at  low  elevations. 

The  present  approach,  which  characterizes  a slant  path  through  the 
atmosphere  by  three  parameters,  which  correspond  to  effective  optical 
depth,  pressure,  and  temperature , is  expected  to  be  more  accurate  than 
the  LOWTRAN  scheme,  which  employs  only  one  parameter,  namely  the 
equivalent  path  length  at  one  atmosphere  pressure. 

In  addition,  the  dashed  line  in  Figures  33-34  shows  transmission 

1 8 

values  derived  by  fitting  Eq.  (71)  to  atmospheric  emission  measurements 
made  from  Mt.  Haleakala  on  the  night  of  January  22,  1976.  In  this  analysis 
it  was  assumed  that  the  atmospheric  absorptance  and  emissivity  were  equal, 
and  that  the  effective  temperature  of  the  atmosphere  was  equal  to  the  local 
ambient  temperature. 

The  slopes  of  the  curves  in  Figure  33  agree  fairly  well.  However 
for  the  8.  3 - 13.  4 pm  window  (Figure  34)  the  slope  of  the  curve  based  on 
emission  measurements  is  seen  to  be  much  less  than  for  those  based  on 
molecular  absorption  calculations.  We  have  not  explained  this  discrepancy. 

To  obtain  emissivity  values  from  the  atmospheric  emission  measure- 
ments an  effective  atmospheric  temperature  must  be  employed.  The  depend- 
ence of  the  derived  emissivity  upon  the  temperature  used  is  presented  in 
Figures  35-36.  Where  data  points  are  missing,  an  emissivity  greater  than 
unity  would  have  been  required  to  account  for  the  measured  intensities. 

For  an  elevation  angle  of  90°  the  spectral  sky  radiance  was  measured 
to  be  1.  6 x 10  watts/cm^- sr- cm  * in  channel  5 and  7.  8 x 10"^  watts/cm^- 
sr-cm  1 in  channel  6.  We  have  already  presented  predictions  of  sky 
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Figure  36.  Emissivity  of  the  Sky  Measured  by  H.  Kent  at  AMOS,  for 
Various  Assumed  Atmospheric  Temperatures 
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radiance  in  channel  6 due  to  H20  lines  for  real  atmospheric  conditions 
in  Figure  19.  The  maximum  predicted  value  is  about  4.5  x 10‘6  watts/cm2  - 
sr-ctn"  If  we  add  the  influence  of  the  H20  continuum  and  C02  lines, 
we  obtain  a maximum  of  about  6.  8 x 10  watts/ cm  - sr-  cm  , which  still 
lies  below  the  measured  data.  The  discrepancy  is  probably  to  be  explained 
by  additional  sources  of  emission,  such  as  cirrus  clouds.  A thin  layer  of 
cirrus  could  easily  have  an  emissivity  ^ of  0.2,  or  a spectral  radiance  at 
230°K  of  1.4  x 10"  watts/cm  - sr-cm  , which  is  sufficient  to  reconcile 
our  calculations  with  the  field  measurements. 

We  conclude  that  real  time  measurements  of  sky  emission  have 
potential  utility  in  determination  of  atmospheric  transmission.  The  present 
study  provides  assistance  in  this  procedure,  especially  in  providing  pre- 
dictions of  the  effective  temperature  of  the  molecular  emitters.  However, 
we  believe  clouds  may  have  a significant  effect,  which  deserves  further 
investigation.  This  study  should  include  the  effect  of  scattering,  which 
reduces  the  transmission  along  a line  of  sight  without  a corresponding 
increase  in  emissivity. 

* * * 

The  basic  conclusion  of  this  work  is  that  direct  or  indirect  determina- 
tion of  the  total  amount  of  water  vapor  overhead  is  a prerequisite  for  the 
quantitative  interpretation  of  infrared  measurements  through  tropical 
atmospheres,  even  from  an  altitude  of  3 km.  The  permissible  elapsed 
time  between  the  primary  and  secondary  measurements  is  estimated  as 
1-2  hours;  a more  precise  determination  of  the  tolerance  would  be  desirable. 
Records  for  10  different  nights  at  AMOS  show  the  local  absolute  humidity 
changed  by  a factor  of  2 or  3 within  20  to  40  minutes  on  three  occasions. 
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